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^ The Fermi Large Area Telescope (LAT) has recently reported the detection of pulsed 7-rays 

from 6 young pulsars (J0631+1036, J0659+1414, J0742-2822, J1420-6048, J1509-5850, and 
^ ^ J1718— 3825), all exhibiting single-peaked pulse profiles [19]. High-quality radio polarization 

'nI" ' data are also available for 5 of these pulsars, allowing derivation of constraints on their viewing 

' geometries. We obtain independent constraints on the viewing geometries of these pulsars by 

' using a geometric pulsar emission code to model the Fermi LAT and radio light curves. We find 

■ fits for the magnetic incUnation and observer angles a and ^ with typical errors of ~ 5°. Our 

results are generally consistent with those obtained by [19], although we do find differences in 
some cases. Our model may lastly provide a framework to constrain the radio emission altitude. 
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1. Introduction 

The Large Area Telescope (LAT) is the primary instrument onboard the Fermi Gamma-ray 
Space Telescope which was launched on 11 June 2008. Fermi LAT is an imaging, wide-field-of- 
view, high-energy /-ray telescope covering the energy range from 20 MeV to >300 GeV. Fermi is 
currently performing an all-sky survey, taking advantage of its large field-of-view and high sensitiv- 
ity to scan the entire sky every three hours [3]. One of the important products to date of this all-sky 
survey is the first Fermi LAT catalogue of 7-ray pulsars [1] produced using the first six months of 
Fermi LAT data, and increasing the number of known 7-ray pulsars from at least 6 [16] to 46. The 
number of 7-ray pulsars continues to increase and is approaching 100 [1 1]. This dramatic increase 
offers the first opportunity to study a sizable population of these high-energy objects. 

We study 6 pulsars (J0631-I-1036, J0659-I-1414, J0742-2822, J1420-6048, J1509-5850, and 
J1718— 3825) detected by Fermi, all exhibiting single-peaked pulse profiles at energies > 0.1 GeV 
within current statistics. These are also detected in the radio band. Both of these properties aid 
in constraining the possible solution space for the respective pulsar geometries. A study of the 
pulsars' geometric parameters, using the rotating vector model (RVM, [14]) for radio polarization, 
as well as a model describing the phase difference between the polarization angle and pulsar Ught 
curve centres [4], has been performed [19]. They constrained the solution space for the geometries 
of the individual pulsars using RVM fits to the radio polarization data, as well as predictions for 
the value of the half opening angle p of the radio beam derived from the radio pulse width (e.g., 
[9]). The aim of this study is to obtain similar constraints on the magnetic inclination and observer 
angles {a and C,) of these pulsars using an independent, multiwavelength approach. We use a 
geometric pulsar emission code (e.g., [17]) to model the Fermi LAT and radio light curves, and fit 
the predicted radio and y-ray profiles to the data concurrently, thereby significantly constraining a 
and ^. This also allows us to test the consistency of the various approaches used to infer the pulsar 
geometry. 

2. JModel 

We use an idealized picture of the pulsar system, wherein the magnetic field has a retarded 
dipole structure [6] and the emission originates in regions of the magnetosphere (refered to as 
'gaps') where the local charge density is sufficiently lower than the Goldreich- Julian charge density 
[10]. These gaps facilitate particle acceleration and radiation. We assume that there are constant- 
emissivity emission layers embedded within the gaps in the pulsar's co-rotating frame. The location 
of these emission layers determines the shape of the light curves, and there are multiple models for 
the geometry of the magnetosphere describing different possible gap configurations. 

We included two models for the 7-ray emission regions in this study, namely the Outer Gap 
(OG) and Two-Pole Caustic (TPC) models. In the both the OG and TPC models (the Slot Gap 
model [2] may serve as physical basis for the latter), emission originates from narrow gaps along 
the last open field lines. In the OG model [5], emission originates above the null charge surface 
and interior to the last open field Une, while the TPC gap starts at the stellar surface. We assume an 
OG emission layer position spanning rove = 0.90 to 0.95, while we set rove = 0.95 to 1 .00 for the 
TPC case (with rove the open volume coordinate which is similar to the polar angle normalized 
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with respect to the polar cap angle [8]). The special relativistic effects of aberration and time-of- 
flight delays, which become important in regions far from the stellar surface (especially near the 
light cylinder), together with the curvature of the magnetic field lines, cause the radiation to form 
caustics (accumulated emission in narrow phase bands) which are detected as peaks in the observed 
light curve [7, 13]. 

We used the empirical radio cone model of [15], where the beam diameter, width, and altitude 
are functions of the pulsar period P, its time derivative P, and the frequency of observation v, in 
order to obtain predictions for the radio Ught curves. (This is different from the approach taken 
by [19], as we have a different prescription for the radio emission altitude and cone half opening 
angle, and we do not use the polarization data.) 

3. Results 

Figure 1 indicates model light curve fits to each of the pulsars, for both the OG and TPC 
models. Figure 2 indicates significance contours in a — j8 space from [19] resulting from radio 
polarization modelling, as well as constraints obtained from the current study. We infer values for 
a and ^ for each pulsar with typical errors of ~ 5° (see Table 1). The latter have been estimated 
by comparing predicted light curves for different [a, 11,) combinations for a 1° resolution, and 
constraining the optimal solution space by eye. The tabulated values are the average values of 
a, C, and /3 implied by the solution contours, while the errors are chosen conservatively so as to 
include the full (non-rectangular) contour We note that once a solution is found at a particular a 
and ^, it is worthwhile to study the light curve at the position a' = and l^' = a (i.e., when the 
angles are interchanged). Such a complementary solution may provide a good fit in some cases (see 
the alternative solutions Usted in Table 1), due to the symmetry of the model (i.e., constant beam 
emissivity assumed) as well the symmetry in pulsar geometry (i.e., |j3 1 remaining almost constant 
under this transformation). 
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Table 1: Best-fit results for the modeled pulsar light curves (all angles are in degrees). 
Alternative best-fit solution (upon interchange of a and C,). 
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Figure 1: Model light curve fits to the 7-ray (panels a, c, e, g, i, and k) and radio (panels b, d, f, h, j, and 1) 
data for each of the pulsars. Blue lines indicate OG fits, red lines TPC ones; dashed lines indicate alternative 
fits (see Table 1). The radio curves result from a conal radio model, but using the same geometry as the 7-ray 
curves. The histograms are Fermi data, while the black lines are 1.4 GHz radio data (reproduced from [19]). 



4. Conclusions 



Generally, we find consistent results compared to tliose of [19], and our mucli smaller solution 
contours usually overlap with their best-fit contours, although we do find different geometeries 
for some pulsars. For better comparison, it should be noted that our model has the symmetric 
property (due to the assumption that the emission is symmetric for both magnetic pole) that when 
one transforms either a' = 180° — a or = 180° — the light curve remains identical, but only 
shifts 180° in phase, as one is then viewing the opposite magnetic pole. This sliift may be neglected 
due to the fact that the light curves are aligned with the radio zero as defined by the data (radio and 
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Figure 2: Independent model constraints on a and ji = ^ — a for each of the pulsars. Black contours are 
reproduced from [19], the green cross signifies their optimal solution. The blue areas are for the OG model, 
and red for the TPC model (solution contours from this work). Contours extending beyond the /3-axis 
indicate that our solutions are not contained in the plotted area. 



/-ray curves are translated in unison so as to preserve the predicted radio-to-/ lag). When doing 
both of these transformations, there is no shift in the light curve, but jS' becomes — j3. 

It is noteworthy that our geometric modelling is able to constrain a significantly better com- 
pared to the constraining power of the RVM. Furthermore, given the small range of jS compared to 
the large range of a encompassed by the contours of [19], it is reassuring to note that we reproduce 
almost all of the j3 -values within errors. The strength of the constraints we are able to derive on 
a and ^ (or j3 ) derives from the combination of having to reproduce the /-ray profile shape, the 
radio peak multiplicity, as well as the radio-to- 7 phase lag. Furthermore, the fact that these pulsars 
have single peaks implies an even smaller possible solution space than for double-peaked profiles, 
which are more common in these 7-ray models [18]. 

Both the OG and TPC models provide reasonable light curve fits, neither being strongly pre- 
ferred. The OG sometimes predict peaks that are too narrow, double peaks with bridge emission 
(between the two peaks) which is too low to reproduce the data, and no off-peak emission. On 
the other hand, the TPC sometimes overpredicts the background, in some cases predicts an ear- 
Uer peak which is not found in the data (double-peaked profile), and generally predicts relatively 
broader peaks. The latter model also exhibits greater variety of peak shapes, and fills a larger region 
of phase space. 

We find that for each pulsar the impact angle p remains the same between the OG and TPC 
geometries, even if the inferred a and differ in each case. Also, j8 is quite small (<20°). This is 
due to the constraints imposed by the 7-ray and radio pulse shapes (as well as their being visible in 
the first place). It is quite striking that j8 < in most cases, so that the observer samples radiation 
coming from above the magnetic axis (i.e., < a). Interestingly, the OG-predicted a generally 
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seems to be larger than that predicted by the TPC geometry. This results from the fact that the OG 
emission layer lies closer to the magnetic axis than that of the TPC model, so to obtain the correct 
radio-to-7 lag, the a must be a little larger. Lastly, we also find that reproduction of the radio-to-y 
phase lag is sensitive to the value of t^. 

Our approach may provide constraints on the radio emission geometry for the different pulsars. 
For example, a much better OG (TPC) fit may be obtained for PSR J0659+1414 if the radio cone's 
altitude is lowered (increased), as this will decrease (increase) the radio-to-y phase lag due to a 
smaller (larger) aberration effect. (Alternatively, changing the position of the emission layer will 
also lead to changes in the y-ray peak phase.) Furthermore, the observed radio profile width may 
be used to constrain the cone width for a given emission height. 

Future work includes applying a Monte Carlo Markov Chain approach for finding optimal 
Ught curve fits in multivariate phase space [12] for these pulsars, constraining the radio beam 
properties, and comparing the new best-fit solutions (using more Fermi data) obtained from the 
different approaches. 
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